1 



ADP-99-38/T375 



Chiral Corrections to Baryon Masses Calculated within Lattice QCD 
Anthony W. Thomas*, Derek B. Leinweber"^, Kazuo Tsushima-^, and Stewart V. Wright^ 

Department of Physics and Mathematical Physics 

and Special Research Centre for the Subatomic Structure of Matter, 

University of Adelaide, Australia 5005 

Consideration of the analytic properties of pion-induced baryon self energies leads to 
new functional forms for the extrapolation of light baryon masses. These functional forms 
reproduce the leading non-analytic behavior of chiral perturbation theory, the correct 
non- analytic behavior at the A^vr threshold and the appropriate heavy-quark limit. They 
involve only three unknown parameters, which may be obtained by fitting lattice QCD 
data. Recent dynamical fermion results from CP-PACS and UKQCD are extrapolated 
using these new functional forms. We also use these functions to probe the limit of 
applicability of chiral perturbation theory. 

1. Introduction 

Chiral symmetry requires that the nucleon mass has the form 

mAr(m^) = m7v(0) + ami + jSml + Inm^ + . . . , 

for small m^, where m^{0), a, /?, and 7 are functions of the strong coupling constant as{fi). 
Recent work has shown that using physical insights from chiral perturbation theory 
and heavy quark effective theory one can derive new functional forms which describe the 
extrapolation of light baryon masses as functions of the pion mass {171^^). These forms are 
applicable beyond the chiral perturbative regime and have been compared successfully 
with predictions from the Cloudy Bag Model and recent dynamical fermion lattice 
QCD calculations. 



2. Analyticity 

By now it is well established that chiral symmetry is dynamically broken in QCD and 
that the pion is almost a Goldstone boson. It is strongly coupled to baryons and therefore 
plays a significant role in the N and A self energies. In the limit where the baryons are 
heavy, the pion-induced self energies of the and A, to one loop, are given by the 
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Figure 1. One-loop pion induced self energy of the nucleon and the delta. 



processes shown in Fig. |I](a-d). We label these by a^N, cttva, o'atv, and cjaa- Note that 
we have restricted the intermediate baryon states to those most strongly coupled, namely 
the N and A states. Other intermediate states are suppressed by the baryon form factor 
describing the extended nature of baryons. 

The leading non-analytic contribution (LNAC) of these self energy diagrams is associ- 
ated with the infrared behavior of the corresponding integrals - i.e., the behavior as the 
loop momentum — > 0. As a consequence, it should not depend on the details of a high 
momentum cut-off, or form factor. In particular, it is sufficient for studying the LNAC 
to evaluate the self energy integrals using a simple sharp cut-off, u{k) = 9{A — k) as the 
choice of form factor. The explicit forms of the self energy contributions for ajy^, a^A and 
so on are given in Moreover, there is little phenomenological difference between this 
step function and the more natural dipole, provided one can tune the cut-off parameter 
A. The self energies involving transitions of — > A or A ^ are characterized by a 
branch point at = AM. 

2.1. Chiral Limit 

The leading non-analytic (LNA) terms are those which correspond to the lowest order 
non-analytic functions of rriq - i.e., odd powers or logarithms of mT^. By expanding the 
expressions given in [IH, we find that the LNA contributions to the nucleon/delta masses 
are in agreement with the well known results of xPT [^,^. 

Of course, our concern with respect to lattice QCD is not so much the behavior as 
ttItt 0, but the extrapolation from high pion masses to the physical pion mass. In this 
context the branch point at = AM^ is at least as important as the LNA behaviour 
near m^r = 0. 

2.2. Heavy Quark Limit 

Heavy quark effective theory suggests that as m^r oo the quarks become static and 
hadron masses become proportional to the quark mass. In this spirit, corrections are 
expected to be of order l/rrig where rUg is the heavy quark mass. Thus we would expect 
the pion induced self energy to vanish as l/m^ as the pion mass increases. The presence of 
a fixed cut-off A acts to suppress the pion induced self energy for increasing pion masses. 
While some dependence in A is expected, this is a second-order effect and does not 
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alter this qualitative feature. Indeed, in the large m^r limit of the equations, we find that 
they tend to zero at least as fast as l/rn^. 

The agreement with both the chiral limit and expected behaviour in the heavy quark 
limit suggests the following functional form for the extrapolation of the nucleon mass [|1| : 

Mn = ttTv + (3Nml + aNN{m^, A) + (TNA{rn^, A) . (1) 

3. Lattice Data Analysis 

We consider two independent lattice simulations of the and A masses from CP-PACS 
and UKQCD Both of these use improved actions to study baryon masses in full 
QCD with two light flavours. We find that the two data sets are consistent, provided one 
allows the parameters introducing the physical scale to float within systematic errors of 
10%. 

We begin by considering the functional form suggested in Section ^ with the cut-off 
A fixed to the value determined by fitting CBM calculations. This is shown as the solid 
curve in Fig. ^. In order to perform model independent fits (i.e. with A unconstrained), it 
is essential to have lattice simulations at light quark masses approaching ~ 0.1 GeV^. 
This fit is illustrated by the dash-dot curve. 

Common practice in the lattice community to use a polynomial expansion for the mass 
dependence of hadron masses. Motivated by xPT the lowest odd power of allowed is 
ml: 

MN = a + (3ml + ^ml (2) 

The result of such a fit for the N is shown in the dashed curve of Fig. ^. The coefficient 
of the ml term, which is the leading non-analytic term in the quark mass, in the three 
parameter fit is —0.761. This disagrees with the coefficient of —5.60 known from xPT 
(which is correctly incorporated in Eq. (|I]), the solid and dash-dot curves) by almost an 
order of magnitude. This clearly indicates the failings of such a simple fitting procedure. 

4. Summary 

In the quest to connect lattice measurements with the physical regime, we have explored 
the quark mass dependence of the N and A baryon masses using arguments based on 
analyticity and heavy quark limits. We have determined a method to access quark masses 
beyond the regime of chiral perturbation theory. This method reproduces the leading 
non-analytic behavior of xPT and accounts for the internal structure for the baryon 
under investigation. We find that the leading non-analytic term of the chiral expansion 
dominates from the chiral limit up to the branch point at = AM ~ 300 MeV, beyond 
which xPT breaks down. The predictions of the CBM, and two-flavour dynamical-fermion 
lattice QCD results, are succinctly described by the formulae derived in The curvature 
around m^^ = AM, neglected in previous extrapolations of the lattice data, leads to shifts 
in the extrapolated masses of the same order as the departure of lattice estimates from 
experimental measurements. 

Acknowledgments 

This work was supported in part by the Australian Research Council. 



4 




Figure 2. A comparison between phenomenological fitting functions for the mass of the 
nucleon. The dotted curve corresponds to using Eq. (0) with 7 set equal to the value 
known from xPT. The three parameter fit (dashed) corresponds to letting 7 vary as an 
unconstrained fit parameter. The solid and dash-dot curves correspond to our preferred fit 
of the functional form of Eq. with A from the CBM and as a fit parameter respectively. 
The lattice data from are CP-PACS (solid) and UKQCD (open), each with a 5% scale 
change. 
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